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We report the observation of the vacuum Rabi splitting of a single quantum emitter by measuring 
its direct spontaneous emission into free space. We used a semiconductor quantum dot inside a 
photonic crystal nanocavity, in conjunction with an appropriate cavity design and filtering with 
a polarizer and an aperture, enabling the extraction of the inherently-weak emitter’s signal. The 
emitter’s vacuum Rabi spectra exhibit clear differences to those measured by detecting the cavity 
photon leakage. Moreover, we observed an asymmetric vacuum Rabi spectrum induced by interfer¬ 
ence between the emitter and cavity detection channels. Our observations lay the groundwork for 
accessing various cavity quantum electrodynamics phenomena that manifest themselves only in the 
emitter’s direct spontaneous emission. 


Cavity quantum electrodynamics (QED) studies the 
interaction between cavity photons and quantum emit¬ 
ters, such as Rydberg ;1] and neutral atoms super¬ 
conducting qubits Q, nitrogen vacancy centers in dia¬ 
mond Q and semiconductor quantum dots (QDs) dm. 
The emitter and cavity photons interact not only with 
each other at a rate of g , but also independently with 
the free-space vacuum field, leading to irreversible radia¬ 
tion at rates of 7 (for emitter) and k (for cavity). When 
9 > K + 7 an d the pure emitter dephasing is negligible, 
the strong coupling regime is achieved and vacuum Rabi 
splitting (VRS) can be observed in the spectral domain. 

In general, the two radiation channels exhibit different 
spectra, which, in principle, can be separately measured 
as the emitter (Sqd(uj)) and cavity (Sc(w)) VRS spec¬ 
tra, as illustrated in Fig. QJa). Moreover, their interfer¬ 
ence (Sj(uj)) should be expected to affect the measured 
spectral shape. However, to date VRS spectra have been 
measured predominantly by cavity transmission, reflec¬ 
tion and emission spectroscopy (2HZ|, since it is brighter 
and easier to access for most cavity QED systems. This is 
particularly the case for QD-based cavity QED systems, 
since most of these possess high cavity leak rates, such 
that k 7 . 

In atomic cavity QED systems, some spectroscopic 
studies on the emitter channel have been performed [lQ| 
but none of these were performed in the single emitter 
strong coupling regime. In contrast, theoreticians fre¬ 
quently discuss the emitter spectra [llUl 8 j|. the impor¬ 
tance of which has been discussed pertaining to the study 
of several intriguing phenomena, including the quantum- 
classical crossover in single atom lasers [ 12 , LL9|, quan¬ 
tum phase transitions in the Jaynes-Cummings Hubbard 
model 0, cavity induced transparency Hzmn in low 
k systems and quantum delayed-choice experiments 0 . 
Importantly, most of these phenomena manifest them¬ 
selves only in the emitter spectra. The significance of 
the interference between the two leakage channels (that 


has a similarity with Fano interference) has also been 
addressed in the literature 0, 123, 24]. 

In this study, we demonstrate an experimental mea¬ 
surement of VRS by detecting direct spontaneous emis¬ 
sion from an emitter into free space, using a single 
InAs/GaAs QD strongly coupled to a photonic crystal 
nanobeam cavity. We show that the simple combination 
of an appropriate cavity design and filtering with both 
a polarizer and aperture enables the extraction of the 
emitter channel contribution from the free space radia¬ 
tion field, which is otherwise dominated by the cavity 
leakage. This technique allows us to isolate and mea¬ 
sure the emitter spectra, Sqd(w), and to compare it with 
Sc(w), as well as to study the interference part, 5/(a;), 
that provides some asymmetry in the VRS spectra. 

We fabricated the photonic crystal nanobeam cavity 
onto a low density quantum dot wafer grown by molecu¬ 
lar beam epitaxy, using a standard combination of elec¬ 
tron beam lithography, dry and wet etching. The details 
of the sample fabrication process can be found in our 
previous publication 0. Figure QJb) shows a scanning 
electron microscope image of a fabricated device. The 
cavity is formed at the center of the air-bridge photonic 
crystal nanobeam by modulating the air hole pattern¬ 
ing period 0, and is designed to have a moderate Q 
factor of 49,000 with suppressed leakage to the vertical 
direction (~ 2% of the total radiation). This suppres¬ 
sion is achieved by tuning the number of reflecting air 
holes [26] and guiding dominant cavity leakage into the 
side waveguide (see Supplementary Section II [0). The 
cavity field distribution overlaid with the cavity design is 
shown in the inset. For optical measurements, the sample 
was placed inside a temperature controlled optical cryo¬ 
stat and was kept at 3.1 K throughout the measurement. 
We use a continuous wave Ti:Sapphire laser oscillating 
at 860 nm to pump the sample with a fixed power of 
7.3 /iW (except when measuring the spectra shown in 
Fig. 09(d)). The pump creates excitons in the QD that 
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radiate as a dipole and finally decay after exciting the 
electromagnetic field in either the emitter channel (Eqd) 
or cavity channel ( Eq ), as illustrated in Fig. [jja). Note 
that the pump power used is well below the saturation 
pump power of the QD under the resonance with the 
cavity mode. 

We used a micro-photoluminescence (fi- PL) technique 
to address the individual sample, with some filters in the 
detection path, as shown in Fig. (He). The pump laser 
light was focused by an objective lens with a numerical 
aperture of 0.65, which is also used for collecting radia¬ 
tion from the sample. Spectra were obtained by a grat¬ 
ing spectrometer equipped with a CCD detector placed 
after the filters (spectral resolution = 13.5 /ieV). A com¬ 
bination of a half wave plate (HWP) at a angle, O', of 
50.5° and a polarizing beam splitter (PBS) enables the 
rejection of the major cavity farfield, which is polarized 
roughly parallel to the x-axis. The aperture, which has a 
diameter of 0.75 mm and reduces the detection numeri¬ 
cal aperture roughly to 0 . 1 , is used to suppress the minor 
cavity farfield which is polarized orthogonally to the ma¬ 
jor one (see Supplementary Section III jSjj). The effect 
of these filters on the measured spectrum of a detuned 
QD-cavity system is shown in Fig. QJd). Without these 
filters, strong cavity mode emission is seen (black curve) 
but it can be largely reduced by the polarization filtering 
(blue curve). Further reduction is obtained by the aper¬ 
ture (red curve), leading to the near-perfect suppression 
of the cavity emission that is required to successfully ex¬ 
tract Sqd{u) near the emitter-cavity resonance (which 
would otherwise be flooded out by Sc(u)). In the fol¬ 
lowing measurements, we rotate the HWP to compare 
the two spectrum components, Sqd(u) and Sc{u), but 
keep the aperture inserted. 

First, we measured the effect of emitter-cavity detun¬ 
ing on the emission spectra at two different HWP angles, 
as plotted in Fig. [2j by tuning the cavity resonant fre¬ 
quency using a Xe gas adsorption technique. When set¬ 
ting a = 5.5°, we obtain spectra that mostly emphasizes 
the cavity emission contribution, Sc{ c^). This situation 
is very similar to almost all previous /i-PL experiments 
based on QD-based cavity QED systems. In Fig. [2fa), 
the cavity mode-like emission, which forms a diagonal 
line in the colorplot, is visible even under detuned condi¬ 
tions due to off-resonant mode coupling that is prominent 
in many QD-based systems [7]. An anti-crossing of the 
two emission peaks is clearly observed as the cavity is 
tuned into resonance with the QD, demonstrating that 
the system is in the strong coupling regime. In contrast, 
we observed a largely different behavior when measuring 
Sqd(u) by setting a = 50.5° (see Fig. [2](b)) . The diago¬ 
nal cavity-like emission line is not visible under detuned 
conditions and the total emission intensity significantly 
decreases as the system is tuned into resonance. This is 
because k 7 in our system and the system energy de¬ 
cays dominantly through the cavity channel when near 



FIG. 1. (Color online) (a) Schematic illustration of the photon 
leakage from our cavity QED system. The direct QD emis¬ 
sion into free space excites E® D through the 7 channel that 
results in the QD spectrum, Sqd(w). The cavity also sponta¬ 
neously radiates and excites E c from the k channel, resulting 
in the cavity spectrum, Sc(w). Interference between the two 
leakage channels also occurs, resulting in Si(uj). Here, k and 
7 include the leakage into both free space and the adjacent 
waveguide modes, and respectively reflect the total leakage 
rate of the cavity and emitter. Some part of the leaked pho¬ 
tons direct toward the detector, conceptually placed above the 
nanobeam cavity, (b) Scanning electron beam micrograph of 
our device. The inset shows the cavity field distribution over¬ 
laid on the cavity design, (c) Schematic of the experimental 
setup. The sample is addressed by a microscope objective 
lens (OL) and pumped by a continuous wave laser reflected 
by the non-polarized beam splitter (NPBS). The collected 
light passes through a halfwave plate (HWP), polarizing beam 
splitter (PBS) and an aperture, and then sent to a spectrom¬ 
eter. The HWP angle, a , is defined as the angle between the 
HWP fast axis and the x axis (which is parallel to the hor¬ 
izontal axis of our optical table), (d) /iPL spectra showing 
the effect of the polarizer and aperture. Without any filter¬ 
ing, bright cavity peak (C), together with a sharp QD peak 
(QD), can be seen (black curve). The rotation of the HWP to 
a = 50.5° strongly suppresses the cavity contribution in the 
spectrum (blue curve). Further cavity mode suppression can 
be obtained by the insertion of the aperture (red curve). 


the resonance condition. Nevertheless, we still observe 
a clear anti-crossing in this detection geometry. A de¬ 
tailed comparison of the two set of spectra in terms of 
intensities, linewidths and peak positions can be found 
in Supplementary Section IV [27]. 


Now, we discuss the VRS spectra at the resonance. 
Figure[3](a) shows a normalized resonance spectrum when 
a = 5.5°. This spectrum is dominated by the cavity chan¬ 
nel, and clearly shows a vacuum Rabi doublet with a 
splitting of 64 /ieV. The emission spectra exhibits some 
asymmetry in that the higher energy emission peak is 
more intense. The origin of this asymmetry is interfer¬ 
ence between the cavity leakage field and a minor con¬ 
tribution of the emitter field, namely due to the term, 
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FIG. 2. (Color online) (a) Cavity channel spectra under var¬ 
ious detunings, a is set to 5.5°. The spectra contain weak 
contributions from the emitter channel leakage (see text), (b) 
Emitter channel spectra under various detunings, a is set to 
50.5°. 


= Re[Eq D Ec ]• This contribution from the emit¬ 
ter field arises because our QD is elliptically polarized 
and its direct spontaneous emission ( Eqd ) can not be 
rejected by the polarization filtering, resulting in the fact 
that we measure Sc{u) + Si (u) (see Supplementary Sec¬ 
tion I and V 0 ). However, by rotating a to 50.5°, we 
can reject the contribution from the cavity leakage and 
can solely measure Sqd(co’), as plotted in Fig. (3^b). The 
doublet now becomes symmetric and exhibits a deeper 
central dip than in the mainly-cavity spectrum. In ad¬ 
dition, the VRS now shows a slightly wider splitting of 
75 geV with a wider spectral distribution in the exterior 
of the doublet. These features are characteristic for VRS 
spectra of the emitter channel under emitter-driven con¬ 
ditions pim. We note that, since the emitter and cav¬ 
ity photon have a complementarity in this linear strong 
coupling regime, we should expect that the cavity VRS 
spectrum under cavity-driven conditions resembles the 
emitter VRS spectrum under emitter-driven conditions 
(see Supplementary Section VI [27]). 


In these plots, the solid lines show the fitting of the 
data to our theoretical model, which carefully consid¬ 
ers the detection geometry (see Supplementary Section 
V 122|). In the model, we explicitly consider the pure 
dephasing of the QD for a better reproduction of the ex¬ 
perimental results. Most of the parameters needed for 
the calculation can be determined experimentally, such 
as \g\ = 41 /xeV, k = 66 geV (Q ~ 20,000) and 7 = 0.28 
/ieV, leaving just three free parameters (emitter pure de¬ 
phasing rate, cavity-free space coupling phase and degree 
of Eqd-Ec field overlap). The experimentally obtained 
parameters accurately reproduce the size of the vacuum 
Rabi splitting measured through the two detection chan¬ 
nels, using formulas derived by C. Cui and Raymer E3, 
namely 2 y/|g | 2 — (ft 2 + y 2 )/8 (= 67 geV) for the cavity 


channel and 2 \/(\g\ A + 2|g| 2 ft(ft + 7 )/ 4) 1 / 2 — ft 2 /4 (= 75 
geV) for the emitter channel. By setting these free pa¬ 
rameters to a set of reasonable values, such as a pure 
dephasing rate of 3 /ieV, we are able to reproduce the ex¬ 
perimental results (including those discussed later). Note 
that, in the model, we fixed the incoherent pumping rate 
of the QD to be 0.065 geV in order to realize the week 
driving condition as in the actual experiments. 

Figure [3^c) shows a logarithmic plot of unnormalized 
resonance spectra for the two detection channels. The 
peak intensity for the emitter side is roughly 100 times 
weaker than the other. The simulation clearly reveals 
the strong asymmetry in the outer sides of the doublet of 
the mainly-cavity spectrum (red). Experimentally, how¬ 
ever, the relatively strong background emission hindered 
observation of this, and the asymmetry is shown by the 
relative intensities of the doublet peaks. 

In FigOI^d), we show a spectrum taken under a roughly 
three times higher pumping power than that used for the 
rest of experiments. Under this strong pumping condi¬ 
tion, an additional peak between the VRS appears in the 
mainly-cavity spectrum (red curve), as indicated by the 
black arrow. This third peak is often experimentally ob¬ 
served in QD-based cavity QED systems and is believed 
to arise from the QD’s spectral blinking and also off- 
resonant cavity feeding 0. On the other hand, the emit¬ 
ter side (green curve) still exhibits the emission doublet. 
The contrast between the two spectra further supports 
the above scenario for the explanation of the origin of 
the triplet. This result showcases the advantages of our 
measurement technique for discussing pure VRS spectra 
without being bothered by the third peak. 

Next, we apply our technique to study the influence of 
the channel interference. We examine how the interfer¬ 
ence modifies the spectra simply by rotating the HWP 
and mixing the two channel contributions in a controlled 
manner. A schematic illustration in Fig. |4ja) explains 
the details of the experiment. Our elliptically polarized 
QD couples to free space by exciting the E® D and E® D 
fields, while the cavity is assumed to contribute only by 
its coupling to the E^ field. By rotating the HWP, we 
can control the field projection angle, O, defined with 
respect to the x axis. Only x polarized light is trans¬ 
mitted by the PBS, such that the projected field after 
the PBS is a mixture of the emitter and cavity chan¬ 
nel contributions and is expressed as E = E£ cos O + 
E® D cos© + £^ D sin©. The measured spectra is pro¬ 
portional to E*E , which, thereby, contains ©-dependent 
interference terms oc Re[E^*E§ D ], Re[E^*E^ D ]. When 
© = 90° + A© with A© ~ 0 and using the fact 
\E£\ \E® D \ ~ \E$ D l the detected field can be ap¬ 

proximated as \E\ 2 ~ \E® D \ 2 - 2AQRe[E® D *EC], and 
therefore we will observe both the pure emitter spectrum 
and a significant contribution from the interference term, 
Re[E® D *E£]- This term is originally zero, since it arises 
from the orthogonally-polarized fields, but becomes finite 
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FIG. 3. (Color online) (a) Normalized spectra at the reso¬ 
nance measured when mostly emphasizing the cavity leakage 
contribution (a = 5.5°), and (b) when accepting only the 
emitter channel contribution (a = 50.5°). (c) Unnormalized 
spectra corresponding to (a) and (b). For these three pan¬ 
els, the balls are experimental data and the solid lines are 
of numerical simulations, details of which are presented in 
Supplementary Section V [27|. (d) Normalized VRS spectra 
obtained under a three times stronger pumping power than in 
(a)-(c). Red (green) balls are of the mainly-cavity (emitter) 
spectrum. The black arrow indicates the emergence of the 
additional peak for the cavity leakage channel. Slight asym¬ 
metry in the QD spectra is due to unintentional detuning 
induced by the increase of the pump power. 


due to the field projection onto the polarizer. 

Figure|4jb) shows a color plot of normalized PL spectra 
at resonance taken whilst varying a from 5.5° to 95.5°, 
which corresponds to a change of 0 from 0° to 180°. A 
complicated change in the spectra can be seen, especially 
around 2a ~ 101° (0 ~ 90°). Numerical calculations for 
different 0s were also performed, and the results plot¬ 
ted in Fig. SJc). The agreement between the two set of 
spectra is remarkable, and the interference effect is high¬ 
lighted when 2a is set to 95° (A0 = —6°) and 107° 
(A0 = +6°), as shown in Fig. 0Jd). Significantly asym¬ 
metric VRS spectra are clearly observed, and the asym¬ 
metry flips between A0 = ±6°. This is readily explained 
by the fact that the interference spectrum is very asym¬ 
metric (as illustrated in Fig. DJa)) and only the sign of 
its contribution flips when the sign of A0 changes. This 
observation clearly indicates the importance of the chan¬ 
nel interference on VRS spectra. We note that the ob¬ 
served spectra contains information of the phase relation¬ 
ship between the cavity, emitter and the free space. By a 
comparison to our theoretical model, we estimate the rel¬ 
ative phase between the cavity and free space field to be 
roughly ^ 0° (see Supplementary Section V E3). This 
capability to know the cavity/free space relative phase 
will be valuable for deeper understanding of open quan- 


FIG. 4. (Color online) (a) Schematic of the interference ex¬ 
periment. The radiation field from the emitter (E § D , E® D ) 
and cavity ( E%) pass through the HWP and are projected 
onto the PBS with a angle, 0. The PBS only allows the 
transmission of x-polarized light and the resulted field be¬ 
comes a mixture of the emitter and cavity fields and is ex¬ 
pressed as E = E^ cos(0) + E® d cos(0) + E® D sin(0). (b) 
Measured VRS spectra for different HWP angles, a. The 
spectra are mostly modified around 2a ~ 101°, where the 
effect of the channel interference is most prominent. Each 
spectra for different a are normalized to its peak, (c) Cal¬ 
culated VRS spectra corresponding to (b). The spectra are 
calculated for different 0 and are normalized as in (b). (d) 
Two asymmetric VRS spectra highlighting the effect of chan¬ 
nel interference. 2 a is set to 95° (A0 = —6° ,lower panel) 
and 107° (A0 = +6° ,upper panel). Solid lines show of the 
numerical simulation. 

turn system theories, as it is often hard to theoretically 
determine its value due to the nontrivial form of the cou¬ 
pling Hamiltonian (28| . 

In summary, we have measured the VRS spectra of a 
single quantum emitter and observed spectral modifica¬ 
tion due to leakage channel interference. We show that a 
simple combination of an appropriate cavity design, a po¬ 
larizer and an aperture is useful for extracting the emit¬ 
ter’s direct spontaneous emission into free space, as well 
as for controlling the interference between the two detec¬ 
tion channels. Our demonstration will provide a means 
to access various intriguing quantum optics phenomena 
that manifest themselves in only the emitter spectrum, 
such as the incoherently-pumped Mollow triplet, which is 
predicted to occur when single emitter cavity QED sys¬ 
tems start lasing H30. In particular, such phenomena 
in nonlinear strong coupling regimes of cavity QED are of 
interest for further experimental studies, since they could 
be hard to reproduce using many-atom cavity QED sys¬ 
tems GHH. We hope that our findings also stimulate 
discussion on novel cavity QED experiments and theories 
that use the emitter direct spontaneous emission, which, 
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for example, may act as an novel feedback channel in 
quantum control. 
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SUPPLEMENTARY INFORMATION 
I. Quantum dot characterization 

Here, we characterize the QD emission used in the ex¬ 
periments under a far-detuned condition, where the cav¬ 
ity is red shifted by 5.7 meV from the QD emission line. 
For the experiments in this subsection, we removed the 
aperture. First, we performed time resolved PL measure¬ 
ment on the QD using a pulse laser source (pulse width 
~ 1 ps) centered at 860 nm and a single photon counter. 
The result is plotted in Fig. lS5l a). The emission life time 
was evaluated to be 2.35 ns (0.28 yeV) by a single ex¬ 
ponential fit after convolving with the detector response 
function (time resolution = 0.45 ns). Since the QD is far 
detuned from the cavity, this value can be considered as 
the spontaneous emission life time without the Purcell 
effect, which is roughly twice slower than the average of 
those for QDs in unprocessed area (1.2 ns). This weak 
photonic bandgap effect compared to two dimensional 
photonic crystals is one of advantages for the usage of 
the one dimensional photonic crystal nanobeam cavity, 
since it has a larger spontaneous emission rate into free 
space, which makes the emitter channel measurements 
brighter and easier. Next, we discuss the polarization 
properties of the QD dipole. We consider that the emis¬ 
sion of the QD arises from the positively-charged exciton, 
which is the dominant exciton species created in our QD 
wafer & This speculation is also supported by the fact 
that the QD emission line shows a linear pump power de¬ 
pendence and does not show fine structure splitting. Ac¬ 
cording to S. Ohno et al. HI , the QD dipole, d a , of the 
positively-charged exciton with a weak heavy hole-light 
hole mixing emits elliptically-polarized light and has a 
form as follows, 

d a = \d a \(xcos9 a + e ±l<l>QD ysm9 a ), (SI) 

where x and y are unit vectors respectively parallel to 
x and y axes. The two signs of c/)qd correspond to two 
possible states in the positively-charged exciton. Fig¬ 
ure [SUb) shows measured PL intensities of the QD emis¬ 
sion plotted as a function of 2x the HWP angle. The 
data is taken under the same continuous wave pump¬ 
ing condition as the experiments in the main text. The 
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FIG. S5. Quantum dot emission properties, a, Time 
resolved PL for the QD emission line. The emission life time 
was evaluated as 2.35 ns (0.28 yeV). b, Normalized PL in¬ 
tensities plotted as a function of 2x the HWP angle. Balls 
(lines) are experimental (fitting) data. Black (red) curves are 
of QD (cavity), c, Measured intensity correlation function for 
the QD emission line. Circles are experimental data. Blue 
and gray lines are respectively of a fitting curve and that af¬ 
ter deconvolution with the detector response function (time 
resolution = 0.7 ns). 


plot also contains a measured curve for the cavity mode 
emission. The curve for QD are fitted with a func¬ 
tion | cos 6 a cos 2a + e ±l( ^ QD sin 9 a sin 2a | 2 , where a is the 
HWP angle. Through the fit, we obtain 6 a = 42.6° and 
(j) QD = 80.8°. Here, we set the x axis parallel to the cavity 
main axis (2a = 21.9°). Finally, we measured the inten¬ 
sity autocorrelation of the QD emission under a different 
detuning condition (the same as that used in Fig. 1(d)). 
The result is shown in Fig. IS5f ch The observed anti¬ 
bunching with the time origin value of 0.17 demonstrates 
that our cavity QED system is working in the quantum 
regime. 


II. Cavity mode design and characterization 

We used a cavity design based on our previous 
work m si. The structure is schematically shown in 
Fig. IS6l ah and has a base lattice constant (a) of 250 
nm, slab thickness ( d ) of 130 nm, slab width (w) of 340 
nm and an air hole radius (r) of 62.5 nm. The refrac¬ 
tive index of the GaAs (n) is assumed to be 3.4. The 
positions of the nearest neighbor holes from the cavity 
center are shifted as illustrated in the figure. We simu¬ 
lated the near- and far-held of the designed cavity by a 
finite difference time domain algorithm provided by Rsoft 
corporation. In order to suppress the light leakage into 
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FIG. S6. Cavity mode characterization, a, Cavity design. The cavity center is indicated by the dashed line. Positions of 
the nearest neighbor airholes from the cavity center were tuned as follows: ao=0.8a, ai=0.822a and a2=0.889a. Nh is number 
of mirror air holes, b, Calculated cavity Q factors for different Nh s. c, Intensity distribution of the E x farfield in the upper 
hemisphere of radiation space. 6 is the polar angle and (ft is the azimuthal angle in the spherical coordinate, d, The same in c 
but for the E y farfield. e, Phase distribution of the E x farfield, in other words, a plot of arg (E x ). 


the vertical direction (perpendicular to the xy plane), we 
reduce the number of mirror air holes and guide the leaky 
light into the nanobeam waveguide. Figure ISCT bi shows 
the simulated cavity Q-factors as a function of the num¬ 
ber of mirror holes. By setting the hole number to 9, we 
have a high enough Q factor of 49,000 while suppressing 
the vertical-top radiation to only 2 % of the total radi¬ 
ation. In this case, ~ 94 % of light is channeled into 
the lateral waveguide. The cavity also supports a small 
mode volume of only 0.29 (A/n) 3 . We then calculated 
the farfield emission pattern by a near-to-far field con¬ 
version. The results are plotted in Figs. IS6l c) and (d). 
The E x farfield is concentrated in the center, while the 
E y farfield is much weaker and found at larger angles. 
Thus, the E y field will be rejected by the insertion of 
the aperture. We also calculated the phase of the E x 
farfield, as can be found in Fig. IS6l ei. Within the nu¬ 
merical aperture of 0.1, the phase of the field is about - 
85°. This phase corresponds to that which the farfield 
acquires when escaping from the cavity. 


III. Effect of aperture 

Here, we discuss how the aperture in the detection 
path affects the suppression of the cavity leakage con¬ 
tribution reaching to the detector. We measured spec¬ 
tra at various aperture positions controlled by a motor¬ 
ized stage. Continuous wave pumping at 860 nm with 
a power of 3.6 fi W is used here and the HWP angle is 
fixed to 50.5° that maximally suppresses the E x cavity 
field intrusion, while the E y cavity field passes through. 
Representative measured spectra are shown in Fig. lS7f a) 
and (b). At the filter center (Fig. IS7t aW strong QD 
emission is observable while the cavity peak is fully sup¬ 
pressed. Meanwhile, the cavity emission becomes visible 
when the aperture is deviated from the center by a few 
hundred /i meters (Fig. IS7f b)). The small peak next to 
the cavity peak is from another irrelevant QD. Next, we 
compared the integrated intensities of the QD emission 
and that of the cavity mode at various aperture positions, 
as in Figs. m c) and (d). Around the position center, 
QD emission becomes stronger while the cavity emission 
is largely suppressed. The cavity field intensity distribu¬ 
tion that is weak around the center resembles that of the 
minor (E y ) cavity farfield (See Fig. l56t c0). This suggests 
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FIG. S7. Effect of aperture, a, Emission spectrum taken 
at the aperture position center, Y=0 pm, Y=0 pm. The 
aperture center is set to the position used for the experiments 
in the main text, where the maximum ratio of the QD emis¬ 
sion to the cavity emission is achieved, b, Emission spec¬ 
trum taken at an aperture position, X=-204 pm, Y 202 
pm. Weak cavity emission is visible, c, Integrated QD emis¬ 
sion intensities plotted as a function of aperture positions. 
Intensities are expressed by the diameter of the circle. The 
inset corresponds to 1000 cps. XY axes here are parallel to 
the xy axes in Fig. 1(c). Note that the aperture has a round 
shape and the diameter is 750 pm. d, The same in c but for 
the cavity emission. The inset corresponds to 100 cps. 


that the aperture mainly helps in suppressing the E y cav¬ 
ity farfleld that can not be rejected by the polarization 
filtering. 


a = 5.5° a = 50.5° 



Cavity detuning (eV) Cavity detuning (|ieV) 


FIG. S 8 . Fitting to VRS spectra presented in Fig. 

2 . Left panels are data for the mainly-cavity emission taken 
with a HWP angle of 5.5°: (top) integrated intensities of 
the respective peaks, (middle) peak energies, and (bottom) 
linewidths. Right panels are for the emitter spectra taken 
with a HWP angle of 50.5°. Red (black) balls are of lower 
(upper) polariton branch. 


nificant differences between the two sets of data are not 
found in their values at each detuning. However, the 
two intensity curves a show clear difference: the total 
intensities of emitter spectra become much weaker as ap¬ 
proaching to the resonance. This can also be confirmed 
by the colorplot in Fig. 2. 


IV. Detuning dependence of VRS spectra 

Here, we show fitting results on the spectra displayed 
in Fig. 2. Fitting was performed by using two Lorentzian 
curves convolved with a Gaussian peak function (FWHM 
~ 13.5 //eV) that represents the spectrometer response 
function. For the emitter spectra (a = 50.5°), if the cav¬ 
ity peak is not found, we reduce the number of Lorentzian 
peak to one. The results (integrated peak intensities, 
peak energies and linewidths) plotted as a function of 
cavity detunings are summarized in Fig. [551 Through the 
fit, we can deduce some of cavity QED parameters. \g\ = 
41/ieV is determined from the value of vacuum Rabi split¬ 
ting and k = 66/xeV is obtained from the fitting linewidth 
for the cavity-like polariton branch under a large detun¬ 
ing condition (see the linewidth plot for the cavity data 
(a = 50.5°)j. For the peak positions and linewidths, sig- 


V. Theoretical modeling 

In this subsection, we derive a set of equations used 
for calculating the spectra in the main text. We con¬ 
sider a cavity QED system composed of a cavity mode 
(at frequency of u c ) and a two level emitter (u a ). The 
system Hamiltonian of this model can be described as 
follows [32j, 

H s = H 0 + Hj (S2) 

Ho = hhj c a[a c + hhj a (i ] G (S3) 

H x =-[ad a + aU a *} ■ E c . (S4) 

Here, a c is the annihilation operator of the cavity photon, 
a the lowering operator for the two level system, d a the 
atomic dipole moment and E c local cavity field vector at 
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the location of the atom. We quantize the cavity mode 
to be E c = ie c (a c — a£), where e c is a vector parallel 
to the cavity field local polarization at the atom location 

and has a maximum length, ^ 2 e<^V > ^ atom is posi¬ 
tioned in the field maximum. Here V is the cavity mode 
volume and e and eo are the relative and vacuum per¬ 
mittivity, respectively. Using Eq. [Si] and rotating wave 
approximation, the interaction part of the above Hamil¬ 
tonian becomes, 

Hi = ihg(cos /? cos 0 a + sin /3 sin 0 a e ±l(t>QD )a\a + ft.c., 

(S5) 

where g = |d a ||e c |/fi, and f3 is the angle between local 
field vector e c and the x axis. The vacuum Rabi splitting 
induced by this Hamiltonian is 2 g — 2g\ cos (3 cos 0 a + 
sin/3 smO a e ±1, ^ QD \. From the fits to the experimental 
VRS data, we obtaine g = 41 peV. The values of 0 a 
and (j)Qu are also experimentally obtained by fitting to 
Fig. IS5f bh Using these values, we obtain possible val¬ 
ues of /? = -87°, -31°, 93°, 149°. For the reproduction 
of experimental data, /3 needs to be -31°. /3 can also be 
set to an equivalent value of 149° (=-31°+180°), which 
corresponds to the case e c —»• — e c , which does not make 
any difference to the reproduced results. 

Next, we consider the system-free space vacuum field 
interaction by introducing the following Hamiltonians, 

HR = hY,“kA,x h k,x> ( S6 ) 

k, X 


dipole interaction, — [ad a + a^d a *] • P/, the emitter-free 
space dipole coupling constant is expressed as, 


h %,\ = ~ ieik ' rae k,\ e k,\ ' da *’ ( S9 ) 

where r a is the atom location. Meanwhile, the form of 
cavity-free space interaction, in other words the expres¬ 
sion for kfe A , may be nontrivial [28|. For the sake of 
later discussion, we also separate Ef into positive and 
negative frequency parts, 


E f - E t + E f > 

E f =i Y, e k,A,x b k,. 
k, a 

E 7 =-*Y. e kAA,> 

k, a 


ikr 


-ikr 


(S 10 ) 


Following a standard procedure with the Born and 
Markov approximations, we obtain a quantum master 
equation that describes the time evolution of the system 
density operator, p, 

^- = - l -[H s ,p] + Lp + L'p, (SH) 

Lp m 2(2crpcr^ — a^ap — pa^a) (S 12 ) 

+ ^( 2 a c pa\ - a\a c p - pa\a c ), 

with 


H S r = 7 k,x b k,x + hac E k k,x h k,x + h - c 

k, A k, A 

where Hr and Hrr respectively correspond to the free 
space reservoir energy and the system-reservoir interac¬ 
tion energy. The latter Hamiltonian is formed under the 
rotating wave approximation, b is the creation op¬ 
erator for the free space field that is specified with a 
wave vector of k and associated polarization A and has a 
frequency of A . 7 ^ A and k^ A are the coupling con¬ 
stants between the free space field and the cavity photon 
and two level system, respectively. Here, the free space 
reservoir includes the modes in the adjacent waveguide. 
Coupling between the system and the waveguide is as¬ 
signed to corresponding wavevectors and polarizations. 
We quantize the free space mode as, 


E f = i J2 e k,x e kAk, 

k, x 


ik r _ h f -ik r \ 
^ °k, X 6 k 


(S 8 ) 


where is the single photon field amplitude for a given 
quantization volume Vq and can be taken as real, namely 

e k x = \! ZcoVq • e k a a un ^ polarizat:ion vector orthog¬ 
onal to k and e l ^ r expresses the spatial distribution. 
Since the emitter and the free space couple through the 


7 = ( S13 ) 

k, x 

K = y^7I-|K fc A | 2 (5(w fe A -w c ). (S14) 

k, x 


These values should match with the experimentally mea¬ 
sured spontaneous emission rate (7 = 0.28 peV) and cav¬ 
ity decay rate (k = 66 peV). The Lamb shifts of the cavity 
and emitter are assumed to be absorbed in their original 
frequencies. We also add other Liouvillian terms: 

L'p = pa — aa^ p — pcrcr^) 

+j p h(2<J^ (jpcr^ <j — cr^crp — pa^ a). 


These are for introducing an incoherent atom pump (P a ) 
and an atomic phase decay ( 7 ^). For all the simulations, 
we set P a =0.065 peV and 7 p h =3 peV. 

Now we consider the state of field released from the 
system into the free space through the system-reservoir 
interaction 33, 34|. The Heisenberg equation of motion 
of the free space field operator is, 

jb Kx {t) = l[Hs + H R + H S R,b kx ] (S15) 

= ~*( w fc,A 6 fc,A(*) + %x a ( t ) + ~ K k,\ ac ^- 
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We can formally integrate this differential equation to 
obtain an expression for b ^ A (t), 

b k ,x(t) = b k , A(0)e-“^‘ 

t 

dt'<j{t')e- iUk ^ t - t ''> (S16) 

t 

dt , a c {t , )e- iuJk ’^ t - t ''> 

The first term corresponds to thermal photons initially 
existing in the free space reservoir and can be neglected. 
The second and third terms respectively correspond to 
the free space fields excited by spontaneous decay of the 
emitter and cavity. Plugging this back to Eq. EM we 
obtain the expression for the positive part of free space 
electric field, E^ (r, t) at a position r : 




where c is the speed of light and 

= icf\d a \ cos 6 a •) (S21) 

hr)y = ief\d a \e ±l(t>QD smO a . (S22) 

Eq. IS2QI has a simple form: coupling constant x wave¬ 
length (27 Tc/uj a ) over distance (r) x single photon field 

amplitude x density of states (D(u a ) = V ^3 ) x system 
operator. 

Second, we apply a similar deformation for the cavity 
leakage contribution and obtain, 

E f + = ef E / k k A e fc,A fc2 sin edkdOdcj) 

A ^ 
t 

dt'a c (t')e- ibJk ^ t - t ' )+ikrcose . (S23) 



£/(r,() = EJ+(r,f) + £:J + (r,t), (S17) 

i J 0 


k, a 

E- 

k, a 


e 5 + = E ‘kAAk [dt'ajty 

J 0 




where we can separately discuss the atomic (E^ + (r,£)) 

and cavity (E^(r,t)) fields. 

We proceed with the theory for the atomic contri¬ 
bution first. We take the atomic position as origin 
(v a = 0) without loss of generality and use a polar 
coordinate that is related to the Cartesian coordinate 
by k = A: (sin 0 cos </>, sin 0 sin </>, cos </>). We only consider 
the farfield detected at a r that is parallel to z-axis ( 
k • r = kr cos 0). At r, the polarization of the farfield 
will be parallel to {x,y}. Moreover, the single photon 
field amplitude can be approximated to be ~ eA 
This is valid as far as the energy of the free space field 
spreads closely around uj a for any {k , A} mode that 
contributes to E~^ . The summation over k can be re¬ 
placed by f dk 3 and dk 3 = k 2 sin OdkdOdcj). Now, 

the atomic decay contribution reads 

E f 7/1 x e k a^ 2 sin OdkdOdcj) 

A f 

t 

dt , a(t , )e-^ h ^ t - t ' )+ikrcoa6 . (S18) 

Following a conventional process [HI, [34] and by using 
Eq. [S9] and the Weisskopf-Wigner approximation, we 
obtain the expression for the farfield (( kr ) _1 << 1): 

E a f + = -^— ( 4a(t--)d a (S19) 

J 87760 r c z c 

1 27rc 

= -* E 8 ex ^* x _ ( S2 °) 

A =x,y 



In order to be in line with the experiment, we ne¬ 
glect the //-polarization contribution in the farfield and 
assume that the cavity mode only couples to the x- 
polarized free space field, k^ x = x - Moreover, we 
assume that k^ x has no dependence on r and takes 
the form, k^ = kf (6)g ((/>). This is for including the 
effect of inhomogeneous radiation pattern of the cav¬ 
ity mode. If f{6 = 0) = o(l), g(cj))d(j) = o(l) and 

fo ^}p-e lkrcos9 d0 ~ o(^), then E^~ in the farfield re¬ 
gion becomes 

1 27TC 

Ej + = -i^xA c R*-^-e f D(oj c )a(t - r ~), (S24) 


where A c represents the radiation pattern effect and 
A c = 4^ Jq 71 g{(j))d(j). k is the cavity-free space coupling 
constant when assuming homogeneous coupling over any 
k. 

Now we separate the x and y polarized contribution 
in the farfield, namely, E^ + (r,£) = E“+ + Ey+ and 

I£j + (r, t) = E^~ and slightly modify the expressions for 
each component vector as, 


E a + = 


ce f ^7rD(uj a ) 

4:hUJ a X 


y/ycos O a xcr{t ), 


(S25) 


+ = Ce/ s in e a e l±<t,QD ya(t ), (S26) 

y 4nc o a r 


E c + = 


ce f \A c \x/7rD(uj c 

AhuJrV 


\—i\/~Ke l9c )xa c (t ), (S27) 


where we redefined the time by setting t — and used 

the fact that e f \d a \/H = yTsfa and l«l = 


We also introduced the cavity-free space coupling phase, 
6> c , by setting kA* = \k\ \A c \e l0c . In experiments, we mea¬ 
sured the field just around a position r || z and used 
a very small numerical aperture (~ 0.1), in which the 
farfield can be assumed to marginally change. Then, the 
measured field, E~jj, will be expressed as a sum of E 
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and after multiplying the detection area. We fur¬ 

ther introduce factors A' and B' in order to take into 
account the difference of the amount of field passing 


through the aperture: 

E+ = A'E C + + B'E a + , (S28) 

and each polarization component reduces to 

E x + = VA^y cos 6 a cr(t), (S29) 

££+ = Vlv^sin 0 a e i± ^ D a{t), (S30) 

E%+ = VB(-i^e- ie *)a c {t), (S31) 


where A and B are parameters that correspond to pho¬ 
ton energy observable within the detection area in the 
unit time due to one photon leakage, respectively for the 
atomic and cavity contributions. The ratio A/B (=2.85) 
can be determined from fitting to experimental data (Fig. 
3(c)). Referring to Fig. 4, the fields after passing through 
the HWP and PBS reads 

Ep D = E c + cos 0 + E%+ cos 0 + E* + sin 0. (S32) 

0 is the angle between the x axis and the rotated E+ 
field. Then, the detected field intensity, Ep D Ep D , can 


be expressed as 

Ep D Ep D = (Ex + cos ©) x h- c - (S33) 

+ (. E a x + cos 0 + E“+ sin 0) x h.c. (S34) 
+ ( E c +E a x ~ + E c x -E«+) cos 2 0 (S35) 

+ (El+E a ~ +E c ~E a +) sin0cos0(S36) 

The first and second term respectively correspond to the 
pure cavity and atomic radiation. The third and fourth 
term are from the interference between x-polarized cav¬ 
ity and emitter, and between x-polarized cavity and y- 
polarized emitter, respectively. Without the insertion of 
the polarizer, the latter term does not exist, so that we 
can call it polarizer-induced interference term. Accord¬ 
ing to the Wiener-Khinchin theorem, emission spectra 
can be calculated by the Fourier transform of the two 
time correlation function, (Ep D (t)Ep D (t + r)). There¬ 
fore, the spectrum in the steady state (t oo), S(u), 
reads 

s M = Km ^ j dre tUT (E^ D (t)Ep tD (t + r))(S37) 
= S c (u >) + S a (u>) + Sn(tsj) + S’/2(w), (S38) 

where 


5 c (w) 
5„(w) 
Sn (w) 
5/2 (w) 


J dTe luJT {a\{t)a c (t + T)), 


Bk cos 2 0 
2n 

Ay | cos 9 a cos © + e l<t>QD sin 9 a sin ©| 


2n 


\J~AB ^/nyp cos 9 a cos 2 © 


2ir 


J dre lUT {a^ ( t)a(t + t )), 

J dre iUT { e - i ^ +e ‘> (J(t)a c (t + r)) + e^ +9 ^ + r))}, 


(539) 

(540) 

(541) 


V AB yfivyp sin 0 a sin 0 cos 0 


27r 


j dre iuT {e~ i ^ +e ^ D) {(j\t)a c {t + r)> + {a\{t)(j{t + r))}. 

(S42) 


p is introduced to account for the field intensity distribu¬ 
tion overlap between the cavity and emitter contributions 
within the aperture (p = 0 , no overlap, p — 1 , perfect 
overlap). 

We solved the master equation (Eq. ISllj) for the steady 
state by setting A. — o, and then calculated the two time 
correlation functions using the quantum regression the¬ 
orem. By Fourier transforming the two time correlation 
functions, as in EqJHUS! we obtain the spectra af¬ 
ter convolving with the detector response function. The 
spectra is calculated as the average of the results for the 
two signs of (j)QD, because our QD can be assumed to 
randomly choose either sign during the measurements. 

Since we treat the interference, the phase of coupling 
constants between the emitter, cavity and free space, 
namely the arguments of g, k and 7 , significantly affect 


on the result. The emitter couples with the cavity and 
the free space through the well-known dipole interaction 
that naturally incorporates the phase information in the 
coupling constants, g and 7 . However, the cavity-free 
space coupling constant, k, in general has a nontrivial 
form [28| and is difficult to determine its phase, which 
is, then, often left ambiguous in the literature. Our work 
here accesses this issue and can experimentally evalu¬ 
ate the phase of cavity-free space coupling, namely 6 C . 
In this discussion, we can neglect the difference in the 
propagation phase between the two detection channels 
(quantified by the distance from the radiation source to 
the detector), since the QD located near the cavity center 
and both of them radiate from a well localized region (<^C 
1 /im) that is much smaller than the emission wavelength 
(^ 1 /im). When reproducing the measured spectra, we 
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had to take \6 C \ < 30° for /? = —31° and the best simula¬ 
tion was obtained when 0 C ~ 0° {6 C ~ 180° for (3 = 149°). 
Therefore, we took 6 C = 0° for simulating all the spectra 
in this work. According to J. T. Shen and S. Fan (35| . 
6 C = 0° indicates that the cavity holds time reversal sym¬ 
metry and a mirror symmetry in the direction of photon 
leakage. It is worth noting that the value largely deviates 
from the coupling phase estimated by numerically solv¬ 
ing Maxwell’s equation -85°) for our cavity structure, 
as plotted in Fig. IS6t e). Further detailed studies will be 
necessary to verify the actual coupling phase between the 
cavity and free space. It will also be interesting to ver¬ 
ify the coupling phase in other radiation directions not 
studied in this work (|| z). 


VI. Emitter-cavity complementarity in VRS spectra 

When under weak pumping, the emitter and cavity 
photon are in the linear strong coupling regime of cavity 
QED. In this regime, an energy quantum injected into 
either the emitter or cavity will exhibit the same fun¬ 
damental vacuum Rabi oscillation in the time domain. 
Therefore, the VRS spectrum of the emitter under weak 
emitter-driven conditions should coincide with the VRS 
spectrum of the cavity under weak cavity-driven condi¬ 
tions. In other words, the emitter and cavity photon 
hold a complementarity [36j. The main difference in¬ 
duced by the different pumping conditions is a different 
initial condition of the vacuum Rabi oscillation, which is 
thus assumed to be the dominant cause of the difference 
between the emitter and cavity VRS spectra under the 
emitter-driven conditions, as observed in the main text. 

In order to verify this, we performed numerical simu¬ 
lations of VRS spectra for two different pumping condi¬ 
tions. For the emitter-driven condition, we used exactly 
the same pumping condition as in Fig.3(a) and (b) in the 
main text, and the results are replotted in Fig. lS9t] a). For 
the cavity-driven condition, we switched off the emitter 
pumping (P a = 0) and introduced an incoherent cavity 
pumping (P c ) through the following Liouvillian, which 
was added to the master equation (Eq. IS11D : 

L"p = L ( 2 a c pa\ - a\a c p - pa\a c ) (S43) 

p 

+-^-(2alpa c - a c a\p - pa c a\), 

where we set P c =0.065 fieV. The simulated spectra for 
the two detection channels are plotted in Fig. IS9f bh As 
expected, the emitter VRS spectrum exhibits a narrower 
splitting and a shallower central dip when the cavity is 
pumped. The remaining asymmetry in the mainly-cavity 
channel spectrum is due to the effect of the interference 
described in the main text. 


a b 



FIG. S9. Comparisons of VRS spectra under two dif¬ 
ferent pump conditions, a, Normalized VRS spectra under 
emitter-driven conditions, as obtained by the numerical simu¬ 
lations (the same plots in Fig.3(a) and (b) in the main text), 
b, Normalized VRS spectra under cavity-driven conditions. 
The spectra were simulated with exactly the same procedure 
as in a, albeit with only the cavity being pumped. In this 
simulation, we set the emitter pumping rate to zero ( P a = 0 ) 
and introduced an incoherent cavity pumping by coupling a 
finite temperature reservoir to the cavity mode. 
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